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Abstract 

We compute corrections to both the isovector anomalous magnetic moment and the isovector 
electromagnetic current of the nucleon to 0{p^) in the framework of covariant two-flavor Baryon 
Chiral Perturbation Theory. We then apply these corrections to lattice data for the anomalous 
magnetic moment from the LHPC, RBC & UKQCD and QCDSF collaborations. 

PACS numbers: 



I. INTRODUCTION 



One of the major challenges in nuclear physics today is to understand the structure of 
the nucleon arising from QCD dynamics. Major progress is made by the continuous im- 
provement of lattice QCD, see e.g. p]. However, in some cases results obtained from lattice 
simulations do not seem to extrapolate naturally to the experimentally known results. A 
prominent example are the lattice data for the electromagnetic form factors of the nucleon, 
a quantity for which high precision experimental data is available [2H5] . It is rather plausible 
that this mismatch is due to the strong dependence on pion mass but quantitative numerical 
proof is still missing. 

The computation of observables on the lattice in general suffers from a number of systematic 
uncertainties: Both lattice spacing and lattice volume are finite, and simultaneously most 
simulations use quark masses that are much larger than the physical ones. Thus the quality 
of lattice results depends on the control over a threefold extrapolation: the continuum ex- 
trapolation (a — i- 0), the extrapolation to the thermodynamic limit {V — )■ 0) and the chiral 
extrapolation {rriq — > m^^^^). 

Therefore it is important to study all observables as functions of the quark masses or equiv- 
alently, as functions of the squared pseudoscalar mass. It is also desirable to be able to 
make a prediction of how large finite volume effects are for a certain observable. A suit- 
able framework that allows one to predict pseudoscalar mass dependence as well as volume 
dependence for baryonic observables is called covariant Baryon Chiral Perturbation Theory 
(BChPT). It has been shown that BChPT provides good chiral extrapolations [HI [7j while 
controlling finite volume corrections [HHH]. 

In this paper we investigate the size of finite volume effects for nucleon form factor calcu- 
lations. We work out the finite volume corrections to the anomalous magnetic moments of 
the nucleon to 0{p^) in two-flavor BChPT. We then confront these corrections with lattice 
data. Using techniques presented in |12j| we also calculate the finite volume corrections to 
the electromagnetic current J^'^. 
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II. FINITE VOLUME CORRECTIONS TO THE ISOSCALAR AND ISOVECTOR 
ANOMALOUS MAGNETIC MOMENT OF THE NUCLEON 



We employ SU{2)f BChPT as introduced in Ref. [13], [Hj. In this effective field theory, 
the elementary degrees of freedom are the pseudo Goldstone bosons (the pion fields) and the 
nucleon fields. In particular, we make use of the so-called modified infrared regularization 
scheme (IR) as described in P, fT5]. This regularization scheme is a modification of the 
infrared regularization proposed in Ref. 

To calculate the anomalous magnetic moment of the nucleon to 0{p^), the effective La- 



grangian 



off 



(1) 



is needed. The contributions ^^t' , -^jvvr have been taken from Refs. [HI HTj. The Dirac 



,(2) ^(1-3) 

and Pauli form factors of the nucleon Fi(Q^) and F2{Q'^) are defined as follows: 



{N{p\s')\J^\N{p,s))=u{p',s') 



2M, 



N 



u{p, s) X T]^!-)], (2) 



where JT"^ = |'U7'^m — ^d'-j^d denotes the electromagnetic current, p, p' are the initial and 
final nucleon momenta, s, s' are the spin vectors and the momentum transfer q is defined 
as q = p' — p. Furthermore, = —q^ denotes the virtuality of the photon and is the 
nucleon mass (M^r = 939 MeV). Here, ?7 is a two-component vector describing the isospin 
content, i.e. t] = (1,0)"^ represents a proton and rj = (0,1)"^ represents a neutron. We 
have Fi{0) = 1 for the proton because JT" is a conserved quantity and -^2(0) measures the 
anomalous magnetic moment in nuclear magnetons. The nucleon form factors are related 
to the isoscalar and isovector form factors via 



FmQ')=F^^{Q'') + F^^>{Q 



-W/n2^ 



F^kQ') = Fi^iQ'')-F^a'iQ 



(p)/n2^ 



-W/n2^ 



(3) 



Thus, we find the following relations for the nucleon matrix elements of the electromagnetic 
current: 



{Nip\s')\JC^\Nip,s)) = uip\s') 
{Nip',s')\j;^\Nip,s))=uip',s') 



r 



^F['\Q^)+t. 



2M, 



N 



l'F[^\Q-)+^a^^^F^\Q^) 



2M, 



N 



u{p,s) X ri^-r], (4) 
u{p, s) X r]^ —r]. (5) 
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(a) (b) (c) 



(d) (e) (f) 

FIG. 1: Feynman diagrams contributing to the anomalous magnetic moments and Ks to 0{p^) 
in BChPT. The soHd hne represents a nucleon, the dashed hne denotes a pion and the wiggly hne 
denotes a photon. Ah pion-nucleon and pion-nucleon-photon vertices stem from the leading order 
Lagrangian, where the pion-photon vertex is contained in ^-^n . 

Here, JT^^^ and J'^-^ denote the isoscalar and isovector currents. The isovector and isoscalar 
form factors are Unear combinations of the quark form factors: 

FSiQ') = FUQ") - FtM) ^ FiAQ")^ (6) 
F[t{Q') = Fl,{Q') + Fl,{Q') ^ Firm. (7) 

Furthermore, we define k„ = F2^\o) = 3.71 at physical m.^. For the isovector form factor, 
the disconnected diagrams, which are usually omitted in lattice simulations, cancel in the 
isospin limit, while for the isoscalar form factor, omitting them adds another source of 
uncertainty. 

In this work, we calculate the finite volume corrections in two different ways: in this section, 
we assume that Lorentz invariance is still intact and we directly calculate the corrections 



for the anomalous magnetic moments. In sect. |IV| we assume Lorentz invariance is broken 
and thus, we have to calculate the corrections to the matrix elements of the electromagnetic 
current because these matrix elements can no longer be decomposed into two form factors. 
The leading order corrections to the anomalous magnetic moments come from the Feynman 
diagrams shown in Fig. [T] and the finite volume corrections are calculated from the same 
Feynman diagrams. We do not impose any restrictions on the time direction, whereas for the 
spatial directions, we assume periodic boundary conditions. If the spatial volume is not 
too small, the finite volume effects stem primarily from pions propagating around the spatial 



box. This leads to the so-called p-expansion, which is valid for small pion masses and large 
spatial volumes {V = L^) such that m^L ^1. To calculate the finite volume corrections 
from these Feynman diagrams, we make the following replacement when integrating over all 
possible loop momenta k = (A;o,k): 

j^snko.k)^-Y.nko,k)-j (^/(^o,k). (8) 

In the sum, k represents a discrete set of momenta that are allowed for a box with periodic 
boundary conditions and a spatial length L, i.e. k = 2iT{nx,ny,nz)'^ /L, where n^^, Uy and 
riz are integer numbers. Using Lorentz decomposition, all appearing tensor integrals can be 
reduced to scalar integrals (see Appendix |B]). Using the definitions 

(9) 
(10) 

(11) 




- ie] [M2 -{p- kf - ie] ' 

where p represents the nucleon four momentum, the 0{p^) result for the finite volume 
corrections can be written as follows: 

SKf\m^,L) = /tf (L) -Kf (oo) 

/2 XydM^^M^dM^ My ^ (12) 

2/2 d d \ „(o)... 1 



(m^,L) = /tf (L)-43)(oo) 



+-^(2— + + (13) 



iJio(L)-iJoi(^) 



^9 A [A 2d 



3/^ [ \ dm? 

am^ am'' 
Here, qa and /tt should be taken in the chiral limit, but to the order these corrections 
have been calculated, using the physical values is accurate. Note that (other than L) no 




FIG. 2: The finite volume correction (5k^ for different values of itIt^L: 2 (solid), 2.5 (dashed), 3 
(dotted) and 3.5 (dashdotted). We have the physical values for qa and Mq = M^v, which is 
valid to the order we are working. 

new parameters have been introduced. For the form of the finite volume functions see 
Appendix [Cj 

At this point, we would like to comment on the use of the simplified fit ansatz 6K,v{mT^, L) = 
aexp{—mTrL) when trying to determine the infinite volume limit, as has been done in [5]. 



This form is motivated by finite volume meson ChPT, i.e. terms of the form (CI), in the 



limit of large ttIt^L and only taking into account the n = 1 contribution. If correct, one 
would expect to find one single value for the correction along the lines of constant m^L. For 
the full expression for 6kv, one finds, however, that the corrections are not constant along 
the lines of constant rriT^L, see fig. [2] Most likely this behaviour is caused by the appearance 
of the mass scale Mq in BChPT. Thus, we would advise to always use these 0{p^) BChPT 
formulae instead. 



III. COMPARISON WITH RESULTS FROM LATTICE QCD 



Using the input parameters given in Table III we can calculate the infinite volume quan- 



tities from the raw lattice data by first rewriting k^, given in lattice magnetons (short l.m.), 
in terms of nuclear magnetons (short n.m.) via the equation [T8] 



/l#phys 

norm lattice N (lA^ 

'^v — j^^jattice ' \^^) 
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and then subtracting the finite volume corrections from this normahzed quantity: 



oo) = <°™(m,, L) - 5kM^, L). (15) 



In Table III we present the input values for the numerical evaluation of the finite volume 
corrections. Since the corrections are proportional to a factor [qa/ f-^)"^, choosing the physical 
values for and instead of the values in the chiral limit amounts to a mere multiplication 
with a common factor. Choosing the physical values for the pion decay constant and the axial 
coupling constant would result in multiplying the corrections with a factor of 0.99. Thus, 
we use the physical values for both qa and fj^, because they are known more accurately 
and it does not change the corrected values k™'"^ significantly. From Table [T] we can clearly 
conclude that the finite volume corrections are considerably large although m^L ^ 3.8 for 
the smallest pion masses. This set of finite volume corrected isovector anomalous magnetic 
moments can be used as input for the second fit. 

The 0{p^) BChPT formula for the anomalous magnetic moment is presented in Ref. fi5\ 119] 
and is of the following form: 



[c, - 16Mom^eU(A) + 4'^ + 4'^ + «(/)] , (16) 
«^'"'^'«r(3,n;-7Af=)log(|^)-3M= 



g\m M, _ ^3^2^2 ^ g^4i ^^^^^^ ( 



(17) 



2 

ml 



Ag\{ce + 1)M2 - gl{5c,ml + 28M^) log 



(TTl \ 



: (bml - 16M2) arccos ^ 



Note that we distinguish between Mq, which is the nucleon mass in the chiral limit, and 
M = Mlm.^), the running nucleon mass which appears in the 0{p^) part of the whole 
O(p^) result. As input for the running nucleon mass we have chosen two different sets of 
LECs which are presented in Tab. |Tlj For the fitting procedures, we have only used 
values for pion masses smaller than 500 MeV, which we consider to be a reasonable fitting 
window, considering that the upper bound estimates for the applicability of BChPT found 
in the literature range from pion masses of 350 MeV to 600 MeV [201422] . depending on the 
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FIG. 3: Fit results for two different sets of input parameters. Tlie shaded area represents the 95% 
confidence level band. Fig. 3(a) corresponds to the LECs labeled Set I in tab. whereas for fig. 
|3(b)[ Set II was used. 

observables considered. The results of these fits are presented in fig. |3} As one can see, the 
finite volume corrected lattice data does not naturally extrapolate to the known value at 
the physical point. This is most likely due to a lack of information in the region where one 
expects a lot of curvature (notice that in the interval 300 MeV < m.^ < 500 MeV all values 
are around k,^ ~ 2.4 within errorbars). 
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IV. FINITE VOLUME CORRECTIONS TO THE ELECTROMAGNETIC CUR- 
RENT OF THE NUCLEON 



For the calculation of the finite volume corrections to the anomalous magnetic moments 
we have assumed Lorentz invariance to still be intact. But since we are doing all calculations 
in a box of finite spatial volume, Lorentz invariance is explicitly broken and thus, Lorentz 
decomposition of tensorial integrals is no longer possible, see Appendix [Cj A recalculation 
of the Feynman diagrams shown in Fig. [T] in a gauge where the polarization of the photon 
is perpendicular to its four momentum {e ■ q = 0) shows that a decomposition of the matrix 
elements into only two different form factors is no longer possible, in accordance with what 
one would expect for broken Lorentz invariance. 



Using Eq. (Dl), one can calculate the corrections to the nucleon matrix element of 
the electromagnetic current for various momenta p and q = (q'o,q)- All formulae cal- 
culated with BChPT are applicable for small q^ (in Ref. |23j it was suggested that 
< 0.1 GeV^). We choose the lowest possible momentum transfers, i.e. q G 
27r/L{(l, 0, 0)"^, (0, 1,0)^, (0,0, 1)^}. Furthermore, we only consider spin projections along 
the z-axis. We then compare the matrix elements of the corrections to the infinite volume 
matrix elements using the 0{p^) BChPT results calculated in [151 IE]- These take the form 



Fi'\q') = l + q'p^^>+Oiq 



Ft\q') = ^. + q'p'^'+0{q'), 



(19) 



where the 0{p'^) results for the slopes p[^^ and p'^' are of the following form: 



(v) 

Pi 



(v) 
P2 



5el(A) 



7g\M^ + 2{5g\ + l)M^\og 



- l5glmlM^ + g\ml{l^ml - AAM^) log 



+ 



rat 



[I5m^ - lAmlMl + IQM^] arccos 



9l 



967T^f^M^{ml - AMI 



2Mo 

-12AM^ + lOSm^Mo^ - ISm^Mg^ 



+ 6(3m^ - 22mlMl + AAmlM^ - 16M^) log 



Mo 



+ 



9l 



487r2/^Mo^m^(4M2 - mlfl^ 



9m° - MmlMl + 2AQm'lM^ - 216m^M° 



+ 16M« 



arccos 



2Mn 



(20) 



(21) 
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To arrive at < 0.1 GeV^ we choose the spatial length of the box to be L = 3.87 fm 
which is significantly larger than the volumes available in current lattice simulations. For 
the numerical evaluation we truncate the infinite sum at In,- 1 = 3 and we have used the 



low energy constants tabulated in Tab. Ill, which are taken from [15J. We have chosen 
the Dirac representation and have normalized the Dirac spinors to u{p,s)u{p,s) = 1. For 
the numerical evaluation of the current matrix elements where we utilize the 0{p^) BChPT 



results, Eqs. (19) and (Dl ), we have used the physical values for gA, fw instead of the values 
in the chiral limit and have also set Mq = Mf^^^, which is accurate to the chiral order all 
these quantities have been calculated in. This then yields the results collected in tab. |V} 
They show clearly that to extract Fi{q^) and F2{q^) one first has to correct the finite volume 
results. Again, one finds that the finite volume corrections to the matrix element of the 
electromagnetic current can be of 0(10%), even for rather large box lengths L. 



V. CONCLUSION 



The results of our work can be summarized as follows: 

1. We have calculated the finite volume corrections to the anomalous magnetic moments 
K,s,v of the nucleon in the framework of SU{2)f covariant BChPT assuming Lorentz 
invariance is still intact. 

2. We have found these corrections to be of 0(10%) for pion masses smaller than < 
300 MeV and box lengths of L ^ 2.5 fm. 

3. The chiral extrapolation of the finite volume corrected lattice data for k„ still does 
not yield the experimentally known value. This is most likely due to the lack of 
information in the region 138 MeV < m-,,: < 250 MeV. Partly, these discrepancies 
might also originate from simulating two instead of three dynamical quark fiavors. 

4. When working in a finite volume, Lorentz invariance is broken. Therefore, we have 
redone all calculations without making use of Lorentz decomposition when dealing 
with tensorial loop integrals. This leads to an expression for the correction to the 
nucleon matrix elements of the electromagnetic current. 
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5. The uncorrected finite volume matrix elements can no longer be decomposed into the 
Pauli and Dirac form factors (see tab. |V]). These corrections are considerably large 
when choosing = rn^^^ and \q^\ ~ 0.1 GeV^, i.e. L = 3.8 fm. 

6. The problem that the lattice data does not naturally extrapolate to the experimentally 
known value for can not be solved by taking finite volume corrections into account. 
However, having data for pion masses smaller than 250 MeV and then applying these 
finite volume corrections should yield a reasonable result. 
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Appendix A: Tables 



TABLE I: Lattice data for botli in MeV and in terms of normalized physical magnetons 
from [2], [3] and [Ij. Also shown are the 0{p^) finite volume corrected values k^™''. To calculate 
we have truncated the infinite sums appearing in ( [l3| ) at n = 20. 



Ref. 


L[fm] 


[MeV] 


Mn [MeV] 


Kt, [n.m.] 


[n.m.] 


m 


2.7 


329 


1154(7) 


2.29(23) 


2.37(23) 




2.7 


416 


1216(7) 


2.63(29) 


2.66(29) 


m 


2.688 


297 


1109(21) 


2.447(99) 


2.554(99) 




2.688 


335 


1172(21) 


2.518(57) 


2.576(57) 




2.688 


403 


1221(21) 


2.508(51) 


2.543(51) 


m 


2.5 


293 


1107.1(111) 


2.456(162) 


2.602(162) 




2.5 


356 


1154.8(80) 


2.475(72) 


2.554(72) 




2.5 


495 


1288.4(80) 


2.339(47) 


2.359(47) 


u 


2.9 


260 




2.608(115) 


2.728(115) 




1.9 


280 




2.505(311) 


2.889(311) 




2.9 


290 




2.466(80) 


2.552(80) 




2.3 


300 




2.877(184) 


3.060(184) 




1.7 


310 




2.396(491) 


2.818(491) 



TABLE II: Low energy constants used as input parameters for the fitting procedure taken at a 
renormalization scale A = 1 GeV. The two sets parameters C2, C3, C4 are taken from [21] while Mq, 
ci and Ci were provided by [25]. 





Mo [GeV] 


ci [GeV-^] 


C2 [GeV-1] 


C3 [GeV-1] 


C4 [GeV-1] 


eUA) [GeV-3] 


Set I 


0.893 


-0.757 


1.71 


-3.64 


2.13 


2.25 


Set II 


0.891 


-0.825 


2.66 


-4.50 


2.45 


2.22 
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TABLE III: Low energy constants used as input parameters for the numerical evaluations. 



9a 


U [MeV] 


5ei(Mo) [GeV-2] 


Se2(Mo)[GeV-2] 


L2695 


92.4 


-L3 


4.47 



TABLE IV: Low energy constants determined from the fit to the pion mass dependent anomalous 



magnetic moment using the O(p^) BChPT formula (16). All values have been obtained at A = 
1 GeV. The first row shows the values and the statistical errors obtained by the fit routine whereas 
the second row shows the values suggested to achieve a 95% confidence level. 



Input 


C6 


e\,,{\) [GeV-3] 


xVd.o.f. 


Set I 


3.61(5) 


0.27(4) 


0.97 




3.61(11) 


0.27(8) 




Set II 


3.71(5) 


0.22(4) 


1.03 




3.71(11) 


0.22(9) 





Appendix B: Basic Integrals 



For all calculations, we have used the following definition for the standard scalar integrals: 



(Bl) 



{2'kY [m? - /2 _ ^e]™ [M2 - {p - If - itf ' 

Here, fh and M are two generic mass parameters, m denotes the number of meson propaga- 
tors and n represents the number of nucleon propagators. By utilizing Feynman parametriza- 
tion, the integrals iJio, -f^oi and Hn can be calculated analytically. They take the following 
form: 



ifio(0,m2,0) = 2m^ 
Hoi{M^,0,0) = 



1 , /m\ 



1 fM 



+ 



p^ - + 



ln(A^ 



^2 



32p27r2 



In ( m 



p^ + - m 



32p27r2 



In (M^ 



(B2) 
(B3) 

(B4) 



32p27r2 



arccos 



- + M2 
2mM 
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TABLE V: Numerical results for the current matrix element, the matrix element of the corrections 
and the matrix element of the electromagnetic current in finite volume. The isospin part of the 
matrix element has been omitted. 



(AT {p, 1) \J^\N {p, i)) (AT {p, 1) \SJ^\N {p, i)) (AT {p, 1) mL)\N {p, i)) 


/ 

(1,0, Of 

\ 


u. (Dooy / 



-0.604959i 



1 


^ n nsc;c:i if? ^ 

— U.UoOOilD 

0.0380878 
0.0738641 - 0.00741085i 
^ 0.0738641 J 




U.DoUioO 

0.0380878 
0.0738641 - 0.612369Z 
^ 0.0738641 y 




(0,1, Of 


u. t oooy / 
0.604959Z 


V / 




' -0.0855116 ^ 
0.0738641 + 0.00741085i 

0.0380878 
^ 0.0738641 y 




' 0.680185 ^ 
0.0738641 + 0.612369Z 

0.0380878 
^ 0.0738641 y 




(0,0, If 


^0.765697^ 



V ) 




^-0.0855116^ 
0.0738641 
0.0738641 

^ 0.0380878 y 




' 0.680185 ^ 
0.0738641 
0.0738641 

(^0.0380878y 





Here we have introduced the quantity 

which contains both the e pole as well as some numerical constants. 7e is the Euler- 
Mascheroni constant. All tensorial integrals of the form 

_a4-^ I ^ {V'.VT....} 

7 (27r)<^[m2-i2_ie]'"[M2-(p-02-ie]"' ^ ^ 



are calculated using the definitions 



— . (B7) 

{2ttY [m? - IT -{P-IW 
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Thus, the functions h[\''^'^^ and if^2 take the following form: 

Hi'^ = ^ (//lo - Hoi + ip' -M' + m')H'ff^ , (B9) 

h[1^ = ^^^^ {2m'Hu - ip' -M' + m')H[^ - H^,) , (BIO) 

= ((2 - d)Ho, + dip' -M' + m')Hi\^ - 2m''H^) , (BU) 

H^^ = ^ + ip' - + m'Ml^ - Ho,) , (B12) 

H[f = (i/o2 + Hi\^ + ip' -M' + m')Hg^ - 2m'HS>) , (B13) 
H[f = ^^^^ ((2 - d)Ho2 + dH[l^ + dip' -M' + m')HS - 2m'Hil^) . (B14) 

Using the following general properties of the scalar integrals 

that are valid for n 7^ and m 7^ respectively, one finds that 



12 " dM^ ' ^ ^' 



where a = 1,2,3. The same principle can be used to calculate -^21^- When employing in- 
frared regularization [T^ or modified infrared regularization [15] in order to obtain properly 
renormalized, scale independent results it is neccessary to calculate infrared integrals. After 
combining the meson propagator with the nucleon propagators using Feynman parametriza- 
tion, one extends the integration 

"1 



dx ^ dx, (B17) 
Jo 

and thus defines the infrared integrals Imn- As described in Ref. [HI [15], we can now 
decompose 

poo pi poo 

/ dx= dx+ dx, (B18) 
Jo Jo Jl 

which corresponds to the decomposition Imn = -f^mn + Rmn, where Rmn are called regulator 
functions. For the functions containing only meson or nucleon propagators, we define 

^mO = HmO, RmO = 0, lon = 0, -Ron = —Hon- (B19) 
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Additionally, the infrared integral In and thus the regulator function can be calculated 
analytically, yielding 



p2 _ J^p p2 _ p2 _ j^2 ^^2 / ^ 

^11 = 2 ^ + qo 2 2 + ifi 2 2 — 

Jm'^m^ - (m2 + M2 - p2)2 / m2 - p2 _ ^2 

arccos ' 



(B20) 



167r2p2 \^ ^4m2p: 

For p2 = ]Vf2 one finds that this expression reduces to 



rr? rr? m} fm\ ^yl 44^2 / m \ 



Appendix C: Finite Volume Functions 

The finite volume functions HiQ{m, L), Hqi{M, L) and H[i\m, M,p, L) for a box of vol- 
ume have been calculated in Refs. [101 EE] and take the following form: 

rhfi{n) 



H^oim, L) = Y: £!%^^i i^nm^L^ ) , (CI) 

n=l * 



?i=i 

i/{?(m, M,p, L) = V ^ / rfxi^o (v/nL2 (^^(^^^ _ (1 _ 2,)p2) + ^2(1 _ _ (^3) 

The function /i(n) represents the multiplicity for each n = n\ + n\ + n\, i.e. the number 
of different combinations for ni, n2 and that yield the wanted n. In our calculation of 
the finite volume corrections to the magnetic moments we have assumed Lorentz invariance 
when dealing with tensorial integrals. 

For the calculation of the corrections to the electromagnetic current, we have assumed 



Lorentz invariance to be broken such that Eqs. (B7) and (B8) can no longer be used. The 



starting point for the calculation of such finite volume functions is again the scalar function 

7 \ ^' ^^4) 



where A is an arbitrary three momentum and is a generic mass parameter. This function 
can be calculated by following the descriptions given in Ref . [12] . This leads to the following 
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expression: 



1 I ? . jji 2 



r(/3)(47r)3/2 

From this expression, the results for the tensorial integrals 



1 . (C5) 



rA,i 



(2^)7 [(k-A)VM2]' 



L3 E 



k*k^ 



(2^)7 [(k-A)V7W2] 



/3- 



can be derived by simply calculating derivatives with respect to A, and Aj. 
results take the form 



rA,j 



1 



d 



/-A _|_ A*rA 



92 



4(/3-l)(/3-2)9A,aA^- 



(C6) 

(C7) 
Thus, the 

(C8) 
(C9) 



+ 



A 



d 



■ d 
A^ — + 5*^ 



2(/3 - 1) [ dKj dAi 



This result agrees with the result derived in Ref. [27]. 



Appendix D: Corrections to the electromagnetic current 

In this section, the results for the recalculation of the Feynman diagrams that contribute 
to the nucleon matrix element of the isovector electromagnetic current to leading one-loop 
order in BChPT are collected. The result can be written as, 



e,{N{p\s')\5JI^{L)\N{p,s)) = u{p\s' 



^ Amp,(A^,L,g,p) 



u{p,s) X r]^—r], (Dl) 
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where we sum over all Feynman diagrams contributing to 0{p^). Thus, the result is of the 
following form: 



Amp(, 



„2 rl 

9a 



dy 





2 ^1 



H,o{m, L) - yM'l,/2{Ml, yp') + M^.-^^h/^iMj, yp') 

d 

Hio{m, L) - yM^h/2{Ml, yp) + M^^- r/i/2(A^i, yp) 



Amp^ = - - (Amp„ + Amp^) 



M 9a , , , , 

+ f „ „^ j ax j dyy 



d{yPi) ' 



4l3/2{Ml p{x, y)) + ?>Mlh,-2{Ml p(x, y)) 



+ 2pi{x,y) 



dpiix,y) 

2^2 /-l rl 

dx dyy 



hi2{M\,p[x,y)) - — ^hi2{M\,p{x,y)) 



dpi{x,yY 



. . 3^p{x,y)Ii/2{Ml,p{x,y)) 

^Jn Jo Jo I 

- y^^i r^ _J zi2{Ml,p{x, y)) - y{e ■ p)li^, sh/2{Ml, p{x, y)) 



dpi{x,y) 



dpi{x,y)dpj{x,y) 



h/2{Ml,p{x,y)) 



Amp^ = 



2/.^ Jo 



dx 



^Ii/2{Ml,xci)+eaj 



d{x^i)d{x<^) 



dpi{x,yy 



I_i/2{Ml,xci) 



Ampg = g\Ampa 
9\M 



fn Jo 



dy 



d 



y{£ ■ p)h/2{Mi,yp) - ^i-Q(^^y^(^i^ yp) 

y{£ ■p)h/2{Mlq{x,y)) 



+ ^M{m^ - q^) / dx dyy 
Jo Jo 



— e 



l3/2{Miq{x,y)) 



dqi{x,y) 

dx I dyy 



(D2) 
(D3) 



(D4) 



(D5) 



(D6) 



-^h/2iMl, q{x, y)) + 3y(£ • y)h/2{Ml, q{x, y)) 



d d 
- Vi^ ■ P)7i gq,^^ y^ hMMl q(x, y)) - ^{x, y)^i g^,(^^ y^ h/2{Ml q(x, y)) 



dqi{x,y)dqj{x,y) 



h/2{Ml,q{x,y)) 



Ampy = 



f2 



(D7) 



In the results above, the following masses 

Ml = x^M^ + xm^, Ml = m'^ - xxq^ , Ml = y^ {M'^ - xxq^) + ym'^ , (D8) 
Ml = fM^ + ym^ - x{l + x)y^q^, (D9) 
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and momenta 



P{x,y) ^y{p - qx), q{x,y) ^yp- xyq, (DIO) 

have been used. They stem from combining the different propagators via Feynman 
parametrization. The Feynman parameter x has solely been used to combine meson prop- 
agators or nucleon propagators, where the parameter y has been used to combine different 
propagators, e.g. a meson propagator with a nucleon propagator. 
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